Corrections to scaling in systems with thermodynamic constraints 
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The theoretical study of phase transitions in classical 
spin models has played an important role in the under- 
standing of the key features of critical behaviour in real 
systems. It is well known that each model like the classi- 
cal Ising, the XY and the Heisenberg model represents a 
certain universality class with a specific set of critical ex- 
ponents. Apart from the spatial dimensionality and the 
range of the interactions a universality class is character- 
ized by the number of components of the order parame- 
ter. In contrast to the case of ideal spin lattice models, 
real systems are exposed to various kinds of imperfec- 
tions (lattice defects, impurities, etc.) which can have 
a significant influence on the true critical behaviour of 
the system considered. In many practical applications 
imperfections in the real systems result in a certain ther- 
modynamic constraint. In the recent years, of notable 
interest were continuum spin fluid models which, 
in particular, are considered as a first step towards the 
modelling of ferrofluids and adsorption surface phe- 
nomena 1^. Several important results, which have both 
theoretical and experimental interest, were obtained for 
such models. For example, it was found that, because 
of the interplay between spin and translational degrees 
of freedom, the phase diagrams for spin fiuids are much 
more complicate 10 -1^ compared with pure liquids 
and could lead to the magnetic ordered state both in the 
gaseous and liquid phases. By applying an external mag- 
netic field one could also shift significantly the locus of 
a gas-liquid transition p^-p^ and change the dynamic 
properties |p^-pl|; both static and dynamic properties 
in the spin fiuid models show the differences from pure 
fluids and/or spin lattice models. Note that in the case 
of spin liquid models we deal, in fact, with a system with 
thermodynamic constraint, namely, fixing the density of 
particles. Therefore, the general question arises how a 
thermodynamic constraint affects the critical singulari- 
ties of the system considered. The answer to this ques- 
tion has been partly given a long time ago by M.Fisher 



p5[ . In particular, it was shown that, due to the con- 
straint, the main singularities in the constrained system 
are renormalized if the critical exponent of the specific 
heat a is positive or remain the same as in the 'ideal sys- 
tem' (without constraint) when a is negative. This ef- 
fect is known in the literature as 'Fisher renormahzation', 
and for many examples of model and real systems with 
constraints this renormahzation of critical exponents has 
been proven. 

In Ref. ||2^ it was found by Monte Carlo simulations 
that the critical exponents of a three-dimensional Heisen- 
berg fiuid model arc in disagreement with the results 
known previously for the lattice model. Similar findings 
were later obtained for three- and two-dimensional (3d 
and 2d) Ising spin fiuids We note that the main 

difference between the cases of 3d Ising and 3d Heisen- 
berg fluids, as systems with an annealed disorder, is that 
according to the Fisher renormahzation ||2^ the critical 
exponents have to be renormalized in the first case and 
remain the same as in the pure system in the second 
one. However, in both cases systematic deviations from 
the critical exponents predicted by theory were found in 
Monte Carlo calculations. We note also that in Refs. 
|p2| , p3| a weak dependence of universal quantities on a 
density of particles n = N/V was observed. From these 
facts we conclude that the critical behaviour in such flu- 
ids with internal degrees of freedom are strongly influ- 
enced by nonlinear crossover effects which hide the true 
asymptotic critical behaviour leading to effective critical 
exponents in power law fits over the restricted tempera- 
ture regions considered. 

In this paper we will focus on the study of continuous 
spin fiuid models and the main problem to be discussed is 
whether the magnetic transition in a fluid spin model lies 
in the same universality class as the corresponding tran- 
sition in the lattice model and, if so, what is the reason 
for such a strong crossover, observed for these systems in 
Monte Carlo simulations 122-^4]. In this paper we will 
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follow the scheme of thermodynamic argumentation de- 
veloped in the original paper of Fisher p5| more than 
thirty years ago. Due to the general arguments used, we 
believe that our results are of interest also to many other 
systems with thermodynamic constraints. 



I. GENERAL BACKGROUND 

Let us consider the grand canonical ensemble with the 
thermodynamic potential Q — Q[fi,h,T] of a continuum 
spin liquid at fixed volume V. The pressure of the system 
is simply connected with ^[fj,, h,T], namely, P = —fl/V, 
and from the thermodynamics we have the well-known 
expression 



dP = nd/i + mdh + sdT, 



(1) 



where n, m, and s are the particles' density, the magne- 
tization and entrophy per unity of the volume, respec- 
tively. The conjugated thermodynamic fields in (|l|) are 
the chemical potential /i, the external magnetic field h 
and the temperature T, respectively. In this study we 
are interested in the ferromagnetic phase transition of 
second order with m being the order parameter. 

The first (and the most important) assumption, which 
has been made by Fisher in Ref. [E5[ and will be also 
used in our treatment, is that the nature of the phase 
transition remains ideal if observed at fixed 'hidden' ther- 
modynamic force which is directly related in our case to 
chemical potential fi. In other words this means that in 
the ensemble with fixed /i the magnetic critical exponents 
of the system considered belong to the same universality 
class as the corresponding lattice model (with the same 
spatial dimensionality and number of components of the 
order parameter and the same properties of the pair in- 
teractions). Of course, the critical temperature Tc of the 
ferromagnetic transition (as well as other nonuniversal 
quantities) will depend on the value of /i, so that in the 
ensemble with fixed /i one has T'c(/i). 

Taking into account this assumption, the thermody- 
namic potential near the phase transition point can be 
written as the sum of two terms 

p = P[T*{T, h, m), /i*(r, h, + AF[T, h, ^i] (2) 
= Psing[r*, h*] + Pr-cs[T*,h*] + AP[r, h,^l], (3) 

where P[T*,h*] describes the critical properties of the 
'ideal' system with temperature T* and magnetic field 
h*- Prcg[T*,h*] and AP[T,h,^] are nonsingular func- 
tions of their arguments with continuous and smooth 
derivatives. The scaled temperature T*{T,h, fi) and the 
scaled field h*{T, h, /z) are also smooth functions, so that 
all the singular features of the system in the vicinity of 
the magnetic phase transition follows from the proper- 
ties of Psing[T*(T,h,fi),h*{T,h,fi)]. Expression (|) can 
be considered as the mathematical formulation of the 
assumption made before. Note that using the ideas of 



scaling theory the form of Pging can be further specified. 
In particular, if the so-called Wegner corrections to 
asymptotic scaling behaviour are neglected, one has (see, 
e.g., ^,^) the following expression 



P- \T* h* 



^* 1 2 — a 



(4) 



where r* = (T* - T*)/T*, z* = h*/\T*f+'f, and f{z) is 
a universal scaling function. 

For the scaled field h*{T, h, /i), taking into account the 
symmetry properties of the magnetic field h we may use 
the following representation 

h*{T,h,fi) ^ h-j{T,h,n) 

with j(T,h,fi) — j{T, —h, ^) > 0. This means that the 
real transition still occurs aX h = 0. 

For the 'ideal' part P[T*,h*] we can write down the 
well known results for the main thermodynamic quanti- 
ties. Namely, in the vicinity of critical point with h* — 
and T* ^ T* one has 
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dP 
df* 



dP 

dm* 
' 'dh* 



h*=0 



So -t- Sir 



A, 



^* 1 1 — a 



ds* 
dT* 



= Co + ClT* -f- Ac 



(5) 



(6) 



(7) 



(8) 



for the magnetization to* , the magnetic susceptibility x* , 
the entropy s*, and the specific heat c* of the 'ideal' sys- 
tem, respectively. The coefficients Am, A^, so, si. As, 
Co, ci, and Ac in (^)-(^ do not depend on r*. The val- 
ues of P, 7, and a are the corresponding critical expo- 
nents of the 'ideal' system, respectively. It can easily be 
proved that the same critical singularities are observed 
at h = h* = for the unconstrained system ^ with the 
fixed chemical potential, when the function T*{T,h, fi) 
has continuous and smooth derivatives. In this case the 
line of critical temperatures Tc{ij,) can be found from the 
equation T*(Tc(/x), 0, /i) = T*, and using the properties 
of the function T*(T, h, /i) it can be shown that the criti- 
cal exponents in the grand canonical ensemble under the 
assumption made above are the same as for the 'ideal' 
system. 

Let us now consider the critical properties of the con- 
strained system with the fixed density of particles n = 
n =const. In this case, if h = h* = 0, one gets from (0): 



dP\ _ ^ fdT* 

f ^[P,,g + AP[T,h,^i]] 
\ 5m 



(9) 



TJi=0 
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The critical temperature in the constrained system 
can be found from the condition = Tc{^c) = Tc{n): 
where the chemical potential fic satisfies the equation 
n{iic, Tc{fJ-c)) = n with the function Tc{^) defined above. 
In particular, this gives the equality T* {T^ ,0, ^c) = T*, 
so that in the vicinity of the critical point [T^, /ic, /i — 0], 
being of our interest, taking into account the properties 
of T*(T,/i,^), one finds 

r*(r,o,/i)^r; + t^A^ + iTAT, (lo) 

with the coefficients and It, where A/z = fi — fj,c and 
AT = T ~ are assumed to be small enough for re- 
stricting of our consideration by the linear terms in (10). 
This expression can be rewritten as follows 



(11) 



with = t^/T* and It = tx/T*. In a similar way, 
taking into account @), one obtains from Eq. (||) the 
following expression for An: 

An = n - n = nir* + ^Jr*!^"" + n2T + n^A/i, (12) 

where r = (T - T^)/T°. The coefficients ni and A,, 
can be expressed in terms of the coefficients si. As, and 
i^, introduced in (R) and (|l0|). For An — 0, combin- 
ing Eqs. (10) and (]l2|), we get the relation between the 
reduced temperature scales r* and r in the 'ideal' and 
constrained systems, respectively: 



with 



and 



T^boT* +blT*\T* 



(13) 



in the constrained system. The critical properties of the 
specific heat follows immediately from (Q) when we use 
Eq. ( |l^ ) for the scaled reduced temperature r*. Taking 
the derivative in (|^) with respect to T, one gets 



ds* 




On the other side, if a > 0, from the expression (^3|), one 
has 

Combining these two expressions, it is easily to show that 
the specific heat Cy remains finite at the critical point in 
the constrained system, and a cusp-like singularity with 
the exponent —a/{l — a) can be found when a > 0. 
Hence, all the main results of Fisher consideration |2^] 
are already reproduced for our specific case. Note also 
that the correction exponents A^, describing the Weg- 
ner corrections to scaling in the 'ideal' system, should 
be renormalized to Ai/(1 — a) if a is positive. In the 
opposite case, when a is negative, the critical exponents 
as well as the correction exponents A^ in the constrained 
system remain the same as in the 'ideal' model. These re- 
sults are known since 1968, when the Fisher's paper p^ ] 
was published. 



II. CORRECTIONS TO SCALING IN THE 
CONSTRAINED SYSTEM 

Let us consider some other consequences which follows 
from the thermodynamic analysis given in the previous 
section. Equations (|l3| ) and (14) may be solved with re- 
spect to T* by iterations, so that we obtain in first order 



— n 

bi ^ An{ <TT^ - n2 



Similar solutions can be found from (|ll|) and ( |l2|) for an- 
other cases, being of interest in experimental situations. 
One may be the case when in the canonical ensemble the 
temperature T is fixed at the critical value T^, r = 0, 
but the density An does change. In this case we obtain 



An = cqt* + cit*\t* 



(14) 



with Co — ni + Ufj^/tf^ and ci = An- The expressions 
( |l3| ) and (O) contain already the central result of Fisher 
analysis ]25| ]: it is seen in (|l^) and ( p^ that, depending 
on the sign of the exponent a, either the first (a < 0) 
or secondfa > 0) term on the right hand side of ( p^ ) 
and/or (|l4| ) is dominant. This directly leads to Fisher 
renormalization of the critical exponents /3 and 7 when 
a is positive. In this case the renormalized exponents 
/3/(l — a) and 7/(1 — a) describe the critical singular- 
ities of magnetization m and magnetic susceptibility x 



T* = sign(r) B|r <; 1 + 61r 1^° + 0(|r|2^° ) ) (15) 



and 



r* =sign(An) C|An|^° |l + c|An|^" -f OdAn]^'^") j, 

(16) 

where the exponents Xq. and A^ depend on the sign of 
the critical exponent of specific heat a, namely. 



T - / , a >0, 

{ 1 , a<0. 



and 



Aq. = < 1 - a 

I —a , a < 0, 



(17) 



(18) 
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respectively. The coefBcients B, b and C, c can be ex- 
pressed by the initial parameters bo, bi and Cq, ci in 
Eqs. ( p^ and (|lj), respectively. For example, one gets 
the following equations for the coefficients B and b: 



B 



,-1/(1-q) ^ ^ 

0^ , a > 0, 

b^^ , a < 0, 



and 



-(l-a)-i6o&r 



1/(1- 



a > 0, 
a < 0. 



(19) 



(20) 



It is worth mentioning that the solutions (|T5|) and ( [l^ ) 
could be used only under some special demands which fol- 
lows form thermodynamic stability conditions (see, e.g., 
[^,|o|). In particular, it could be shown |29| that for 
a > a second order phase transition in the constrained 
system is observed only if bi in Eq. ( |l3|) is positive. For 
&i < in the vicinity of phase transition (small r*) the 
constrained system becomes unstable and the second or- 
der phase transition transforms to a first order transi- 
tion. Note that the condition 6i = gives the position 
of a special tricritical point with the critical exponents of 
the 'ideal' system. Similar consideration may be applied 
for the opposite case with a < 0. Therefore, let us as- 
sume that the parameters bo and bi (cq and ci) are such 
that a second order phase transition still exists in the 
constrained system, and the main question then arises 
what are the singular properties of this transition in the 
nonasymptotic region. 

The first important conclusion can be drawn from 
Eqs. (|i|) and (0). It is seen in Eq. ( p^ ) that, indepen- 
dently on either the Fisher renormalization changes the 
critical exponents in the constrained system or not, the 
new corrections to scaling with the exponent appear 
due to the constraint. The exponent of these corrections 
are proportional to \a\ and, because this value is smaller 
for the Ising, the XY and the Heisenberg universality 
classes compared to the Wegner corrections A^, one can 
expect significant contributions from the new corrections 
to scaling just in the pre- asymptotic region. For exam- 
ple, the magnetization m in the constrained system when 
a < can be written as follows: 



M ^ A„,\Tf { I + ai\T 



+ 0(|t|2|"I) I, (21) 



where Ai = uji/ is well-known Wegner correction (see, 
e.g., iQ). Because of usually \a\ < Ai, one can con- 
clude, therefore, that Fisher correction with the exponent 
|a| is dominant in (^l|). In a similar way the case a > 
may be considered. It is evident for both cases that the 
width of asymptotic region in temperature scale, where 
we can restrict the description of critical singularities by 
considering the main critical exponents only, is reduced 
significantly. 

The new corrections to scaling (see (|l^) and (|l^)) with 
the exponents (p^) have appeared in our treatment as the 



result of the temperature rescaling due to the constraint. 
Similar idea was recently used by Krech p9| for the study 
of spin lattice models with constraints within the finite 
size scaling technique. However, we note that the tem- 
perature rescaling is not the only source for corrections of 
that type. For proving of this statement let us consider 
in more details the thermodynamic transformations from 
the ensemble with fixed to the ensemble with fixed n. 
But before proceeding further we derive some additional 
relations to be useful for the subsequent calculations. 
From Eq. (|l|) one derives 



^) -h+n,(l-a)|r*|-"]f^l +n,. 



(22) 



The derivative dr* /dfi can be easily found from Eq. (|1 
and inserting we obtain 



(p) =h + ci(l-a)|r* 



(23) 



where the coefficients cq and ci are the same as in 
Eq. (|l|). RecaUing now that the derivative ( |2^ ) is di- 
rectly related to the compressibility kt of the system con- 
sidered {{dn/d^)rp j^^Q = v?kt), we can conclude that for 
a > a weak divergence appears in the isothermal com- 
pressibility in the both ensembles. The critical exponent 
describing such singular behaviour in the constrained sys- 
tem is equal to the renormalized index of specific heat 
a/ (1 — a), and 
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with the parameters kJ^ and — i^n?ci{l — a)B. This 
shows the main difference to the case of the specific heat, 
discussed above. In the opposite case, when a is nega- 
tive, the compressibility in the both ensembles displays 
the same cusp-like singularity with the exponent a as 
is observed in the specific heat behaviour. Hence, our 
second conclusion is: for a > in the constrained sys- 
tem there are some thermodynamic quantities (e.g., the 
compressibility in our case) which are asymptotic weak 
divergent with the critical exponent a /{I — a). In gen- 
eral, these quantities can be expressed by the second or- 
der derivatives with respect to a 'hidden' thermodynamic 
force. 

The result obtained above is important also for the 
study of the corrections to scaling in thermodynamic 
derivatives. In order to illustrate this statement let us 
use the well-known thermodynamic relation 



XT,, 



/9M\ _ /dM 
\'dh)T,n~\'9h 



fdM_y fdnV 



(24) 
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which gives the connection between the magnetic suscep- 
tibihty in the canonical and the grand canonical ensem- 
bles. For the first term in the right hand side of (|2j) the 
rescaling formula (^) for temperature dependence can be 
used. In the vicinity of the critical point, the contribu- 
tion from the second term in (|24|) can be easily estimated 
by using Eqs. (||), (pi]), and (^3|). For the case of posi- 
tive a this term produces the contributions proportional 
to |t-*|2/3-2+q^ |^*|2/3-2^ ^^^^ Applying the hyperscaling 

relation 2/3 -I- 7 = dv, one obtains for these terms the 
following estimations |r* j"''', Jr* | etc., and, taking 
into account the expression (iq), one can conclude that 
an additional correction term with the exponent a/{l—a) 
appears in the canonical magnetic susceptibility. Similar 
consideration could be used for the case of negative a. 
Hence, the third conclusion is that in addition to the 
rescaling mechanism, given by ^idj), there exists another 
source for the appearance of the same type of corrections 
to scaling (proportional to a), due to the Legendre trans- 
formation. In particular, this means that the new cor- 
rections to scaling can not be included in the standard 
finite size scaling technique just by simple rescaling of 
the reduce temperature (Qq) as it was proposed in p9[ . 

The results obtained can be easily generalized for: (i) 
the case of more complicate constraint, formulated in the 
form F[n,T,^i] = 9 =const (see, e.g., [|9|j30|); and (ii) 
the case, when several 'hidden' thermodynamic forces 
exist in a system. In the case (i), as it was shown in 
[^,0, the generalized constraint will modify the coef- 
ficient 60 and 61 in (^3|) and change the thermodynamic 
stability conditions, so that the phase transition of the 
second order may transform to a first-order transition. 
In fact, the constraint function F[n,T, /z] = =const 
allows to reformulate all the results obtained in some 
'mixed' ensemble with 9 considered as a thermodynamic 
parameter, but the conclusions made above will be still 
valid. Some exceptions may be noted for special choices 
of the constraint function F[n,T,fi\ (e.g., in the trivial 
case F[n,T,fi] = f[T,^l]). 

In case (ii), when several 'hidden' thermodynamic 
forces {fJ.i} with i = 1, 2, . . . , / exist in a system and some 
of them are constrained by external conditions, we can 
applied the scheme described above by performing con- 
secutively Legendre transformations, starting from any 
constraint condition. It is evident that, for a > in 
the 'ideal' unconstrained system, after the first transfor- 
mation we will find renormalized Fisher exponents and 
the new corrections to scaling, so that the critical expo- 
nent of specific heat will change sign {a — > ~a/{l — a)). 
This means that the critical exponents obtained are not 
changed any in further transformations, because of neg- 
ative value of the specific heat exponent — a/(l — a). 
Any next Legendre transformation may change only the 
amplitudes of the singular terms. Note also that the cor- 
rections to scaling, which will appear after the second 
(third and so on) Legendre transformation, will have the 
same exponents (multiple to a /{I — a)) as found after 
the first step. The case with a < in the 'ideal' uncon- 



strained system is even more trivial, and after the first 
transformation we will find the same ('ideal') critical ex- 
ponents with the new corrections to scaling proportional 
to \a\. This picture will not change in further transfor- 
mations. Hence, we conclude that if in the constrained 
system with arbitrary constraints a second order phase 
transition is observed, one has to expect that the critical 
exponents, describing this transition in the asymptotic re- 
gion, are defined by the rules, established by Fisher Mi5 
and the leading corrections to scaling are given by Hit 



III. DISCUSSION AND CONCLUDING 
REMARKS 



We end with a few concluding remarks: 

(i) The results presented are quite general and impor- 
tant for many applications. It is worth to note that con- 
strained equilibrium systems are widely studied in the 
literature, but only a few examples could be mentioned 
when the Fisher corrections to scaling were included into 
consideration (see, e.g., ||3l|-^). In general, such cor- 
rections may affect significantly the critical behaviour 
of compressible magnets [^-33|, systems described by 
Hubbard model [Q, binary and multicomponent mix- 
tures |36l|37|], magnetic spin liquids [2^ [2^, and ''He-^He 
mixtures | |38| , |39t . 

(ii) It seems interesting to study more carefully the 
critical properties of magnets with quenched disorder. 
The motivation is the following. One of the approaches 
to the quenched systems is based on the idea of Morita 
pO|-^. In this approach a quenched system is consid- 
ered as a quasi-equilibrium system with additional forces 
of constraints, keeping the quenched system in its equi- 
librium state and fixing the moments of random distri- 
bution that defines the disordered state. Therefore, the 
quenched system can be studied by means of the equi- 
librium theory, using some results presented in this pa- 
per. In this connection we point out some known results, 
which support this view: (a) the so-called Harris crite- 
rion 1^^, proposed for the quenched systems, is formally 
similar to the Fisher criterion [ p5| , concerning the diver- 
gence of the specific heat in the constrained system; (b) 
in Ref. ||50|] within the e-expansion scheme it was proved 
analytically (using the replica trick) that the leading cor- 
rection to scaling in a randomly diluted inhomogeneous 
0{m) Heisenberg model is a term with the correction ex- 
ponent —a, as it is expected for the constrained Heisen- 
berg model; (c) in Ref. |43| it was shown that if a > in 
the pure system, the traditional renormalization group 
flows, describing the disorder-induced universal critical 
behaviour are unstable with respect to replica-symmetry 
breaking potentials, found in spin glasses. 

(iii) There were several studies of constrained systems, 
performed within the renormalization group technique 
||l|,|l]j42|j3|,||-|3. In general the results obtained sup- 
port the conclusions made within the thermodynamic 
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treatment: Four fixed points (gaussian, spherical, pure 
model and renormalized) were found within the renor- 
malization group theory, and which of these is the stable 
one (pure or renormalized) depends on the sign of the 
exponent a in the 'ideal' system. However, we note that, 
because the stability condition depends on a, it is signifi- 
cant to reconsider these results in higher approximations 
and, in particular, to investigate the effective critical ex- 
ponents for the most typical renormalization group flows. 

(iv) New corrections to scaling, because of smallness 
a for the Ising, the XY, and the Heisenberg models, 
have to be taken into account within standard finite size 
scaling technique In particular, we suppose that 

the methodological problems found for three- and two- 
dimensional Ising model fluid (2^,^J] are directly con- 
nected with the strong crossover effects caused mainly 
by Fisher corrections to scaling. 
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